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Understanding the loss of magnetic order and the microscopic mechanisms involved in laser induced 
magnetization dynamics is one of the most challenging topics in today's magnetism research. WhUe 
scattering between spins, phonons, magnons and electrons have been proposed as sources for dissipation of 
spin angular momentum, ultrafast spin dependent transport of hot electrons has been pointed out as a 
potential candidate to explain ultrafast demagnetization without resorting to any spin dissipation channel. 
Here we use time resolved magneto-optical Kerr measurements to extract the influence of spin dependent 
transport on the demagnetization dynamics taking place in magnetic samples with alternating domains with 
opposite magnetization directions. We unambiguously show that whatever the sample magnetic 
configuration, the demagnetization takes place during the same time, demonstrating that hot electrons spin 
dependent transfer between neighboring domains does not alter the ultrafast magnetization dynamics in 
our systems with perpendicular anisotropy and 140 nm domain sizes. 



Since the first experimental proof of ultrafast demagnetization in 1996 by Beaurepaire et al. \ laser induced 
magnetization dynamics has been the subject of abundant theoretical and experimental studies. The 
motivation is driven by the twofold challenge of understanding magnetization dynamics in a strongly 
out-of-equilibrium system and optically controlling the magnetic state of a nanostructure on the shortest time- 
scale. Apart from the potential applications in magnetic storage media and other spintronic devices, the under- 
standing of this phenomenon remains incomplete and no consensus on the importance of the different 
microscopic mechanisms involved in the ultrafast magnetization dynamics has been achieved so far. 

Indeed, different mechanisms have been proposed to explain the loss of magnetization following a strong laser 
pulse excitation. Initially, they were based on the conservation of angular momentum during the demagnetization 
process. Some of them relied on the coupling between the photon field and the spin bath^'^ while others were 
founded on spin-flip processes occurring when electrons scatter with electrons*, phonons^ or magnons"''^. 
Recently, the role of optically induced ultrafast spin transport on the magnetization dynamics started to be 
addressed. It was first shown by Malinowski et a/."" that hot electrons spin transport enhances and speeds up 
the demagnetization in a spin valve structure as compared to a single magnetic layer. Thereafter, Battiato et al. 
developed a semi-classical model based on optically excited spin polarized electron transport in the superdiffusive 
regime". The spin- dependent lifetimes and velocity of optically excited hot electrons result in spin currents 
within a ferromagnetic layer, inducing a modification of the magnetization profile within the layer on the 
femtosecond timescale. Consequently, it was argued that superdiffusive spin transport may alone be responsible 
for the ultrafast magnetization dynamics'". The existence of laser induced hot electron spin transport was then 
clearly observed by Melnikov et aV. 

The recent advent of techniques combining femtosecond time resolution and element-specific measurement of 
the magnetization such as synchrotron radiation using femtoslicing^ ", free-electron lasers^"'^' and high-harmonic 
generation^^'^^ sources has allowed tremendous progress in this field of research. For instance, Mathias et al. have 
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shown that two entities in a strongly coupled ferromagnetic alloy can 
have their own demagnetization dynamics^**. In the case of ferrimag- 
netic alloys, the different dynamics of both sublattices has been 
shown to be an important ingredient in laser induced magnetization 
reversaP^. Regarding optically excited hot electrons, it was shown 
that hot electron spin transport on the femtosecond timescale could 
not only enhance the demagnetization but could also increase the 
magnetization of a buried layer in a magnetic multilayer depending 
on the relative magnetization orientation in the different layers""". By 
controlling the inter-layer spin current, it was demonstrated that 
both spin transport and spin-flip scattering processes contribute to 
the ultrafast magnetization dynamics^'. However, in the case of single 
magnetic layers, the importance of ultrafast spin transport remains 
unclear^". While all these works dealt with spin transport taking place 
along the multilayers thickness, an interesting experiment was car- 
ried out by Pfau et al?'' and Vodungbo et air'" who studied the 
influence of superdiffusive transport on the demagnetization in the 
presence of nanoscale magnetic domains. In this case, the samples 
are chemically homogeneous and the boundaries between magnetic 
domains are formed by magnetic domain walls (DW). It was shown 
that hot electrons spin transport induced a modification of the mag- 
netic form factor due to spin accumulation around the DW, resulting 
in a change of the magnetic scattering pattern^'. By studying the 
evolution of the demagnetization time (xm) as a function of the 
magnetization quenching, Vodungbo et al. concluded that superdif- 
fusive spin currents were responsible for the constant demagnetiza- 
tion time independently of the laser fluence''". 

In this paper, we report on the influence of the domain configura- 
tion in ferromagnetic structures on the ultrafast magnetization 
dynamics. We show that the demagnetization time increases with 
the laser fluence when the sample magnetization is saturated, simi- 
larly to what is usually observed in transition metals. Furthermore, 
our results reveal that the demagnetization time does not evolve 
when the sample magnetic configuration goes from saturation to a 
multidomain state, questioning the role of hot electrons spin trans- 
port on the ultrafast demagnetization in such a system. 

Results 

Multilayer fabrication and static characterization. The samples 
used in this study consist of a two different multilayers [(Co 
(0.4 nm/Pd (0.6 nm)]x2o (CoPd multilayer) and [(Co (0.5 nm/Pt 
(1.0 nm)]xi5 (CoPt multilayer) grown at room temperature by dc 
magnetron sputtering. The static magnetic properties were analyzed 
by magneto-optical Kerr effect and magnetic force microscopy. 
Figure lb and Ic show hysteresis loops measured on both samples 
with a magnetic field applied perpendicularly to the sample plane. 
From these measurements, it is clear that the CoPd and CoPt 
multilayers exhibit an out-of-plane anisotropy with a domain maze 
structure at low applied magnetic field showing alternating up and 
down magnetization (Fig. la). Indeed, when the field is reduced from 
saturation, a large jump in the Kerr signal is observed even though 
the magnetic field is not reversed. This jump is due to the onset of the 
domain formation. The onset field values are slightly different for the 
CoPd multilayers (80 mT) and the CoPt multilayers (50 mT) due to 
their different saturation magnetization and effective anisotropy. 
Magnetic force microscopy images obtained for both samples are 
depicted in the insets of figure lb and Ic. Both samples show 
similar domain configurations with domain sizes of 145 nm and 
132 nm for the CoPt and CoPd multilayers, respectively. Further- 
more, by modifying the amplitude of the applied magnetic field, it is 
possible to control the magnetic structure of our samples, i.e. the ratio 
between the size of domains with opposite magnetization direction. 

Magnetic domain structure and ultrafast demagnetization. To 

investigate the role of hot electron spin transport between nei- 
ghboring domains, we performed a systematic study of the 




-300 -200 -100 0 100 200 300 

u H (mT) 




Figure 1 | (a) Schematic representation of the sample magnetic domain 
structure. Blue and yellow areas represent magnetic domains with opposite 
magnetization direction perpendicular to the film plane. The blue and red 
curved arrows schematically show the probe and pump laser pulses, 
respectively. Magnetic hysteresis loop showing the variation of the Kerr 
rotation with the applied magnetic field for the (b) CoPt and (c) CoPd 
multilayers. The fiUed circles represent values of the Kerr rotation at a given 
field for which time resolved measurements were performed.(insets) 
Magnetic force microscopy images measured at zero applied field on the 
(b) CoPt and (c) CoPd multilayers. The images size is 5 X 5 )rm^. 

demagnetization occurring in our multilayers after a laser pulse 
excitation for a large panel of fluences and applied magnetic fields. 
The magnetization dynamics was studied by polar time resolved 
magneto-optical Kerr effect (TRMOKE) measurements. The sam- 
ples were pumped by an infrared laser pulse at 800 nm wavelength 
with a pulse duration of 50 ± 5 fs. The probe beam was frequency 
doubled (400 nm wavelength) in order to minimize optical artifacts 
such as dichroic bleaching". In this configuration, the measured Kerr 
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rotation is proportional to the out-of-plane component of the 
magnetization. The measurements were performed for a wide 
range of fluences (0.35-5.6 mj/cm^) and for different amplitudes 
of the magnetic field applied perpendicularly to the sample plane 
(0-300 mT). 

To start with, we performed time resolved experiment as a func- 
tion of the laser fluence under a strong applied magnetic field to 
saturate the samples magnetization. The results are presented in 
fig. 2. Both multilayers show similar demagnetization curves with a 
fast loss of the Kerr signal occurring within the first hundreds of 
femtoseconds and followed by a slow magnetization recovery after 
few picoseconds (Fig. 2a,b). The absorption of a strong laser pulse 
results in a quasi-instantaneous increase of the electronic temper- 
ature. The energy is rapidly transferred to the spin bath leading to the 
fast drop of the magnetic signal. Equilibration between the electrons 
and phonons baths due to the electron-phonon interaction leads to 
the thermalization of the spin bath and a partial recovery of the 
magnetization. It is interesting to note that the time needed to 
demagnetize the sample increases with the demagnetization ampli- 
tude and so, with the fluence. Fig. 2c reports the evolution of the 
demagnetization time obtained by fitting the data using the three 
temperatures model (see Methods) as a function of the demagnet- 
ization amplitude''^^. We can clearly observe a large increase of the 
demagnetization time with the quenching of the Kerr signal for both 
CoPd and CoPt multilayers. This trend is similar to what has been 
previously observed in transition metals thin films'. The existence of 
a difference in demagnetization time for both multilayers might 
come from the difference in relative thickness between the Co and 
the Pd or Pt layers inducing different spin flip probabilities. 

In order to study the influence of the magnetic domains config- 
uration on the ultrafast demagnetization in these multilayers, we 
performed TRMOKE measurements for different applied magnetic 
fields. The measurements were performed at a fixed laser fluence to 
avoid spurious variation of the demagnetization time. Therefore, it is 
possible to carefully control the magnetic domain structure while 
measuring the ultrafast demagnetization. Fig. 3a shows measure- 
ments obtained on the CoPd and CoPt multilayers and for different 
values of the applied magnetic field. While above 80 mT (50 mT) the 
CoPd (CoPt) multilayer is completely saturated, reducing the applied 
magnetic field leads to the appearence of a magnetic domain network 
configuration. Time resolved measurements performed for these 
different applied fields clearly look identical, showing no influence 
of the domain structure on the ultrafast demagnetization (Fig. 3a). To 
get a quantitative estimate of the role played by hot electrons spin 
transport in our systems, we extracted the demagnetization time 
from these measurements. The evolution of the demagnetization 
time as a function of the magnetic configuration is plotted in 
Fig. 3b for CoPd and CoPt multilayers, for a laser fluence of 
4.9 mj/cm^. We used the largest fluence available to maximise the 
number of excited electrons and because a higher fluence results in a 
longer demagnetization time as shown in Figure 2c. The CoPt and 
CoPd demagnetization times, i.e. T^°^' = 200fs and T^°^'' = 180fs 
remain constant independently of the magnetic domain structure, 
indicating no or small contribution of the hot electrons spin trans- 
port between neighboring domains. 

Discussion 

To understand the role played by the hot electrons displacement 
between neighboring domains on the ultrafast demagnetization 
dynamics, we need to take into account that hot electrons transport 
properties are spin dependent. Indeed, in a magnetic material, spin- 
minority (minority) and spin-majority (majority) hot electrons pos- 
sess different lifetimes and velocities'"'". Generally, minority hot 
electrons have a shorter lifetime and a lower velocity compared to 
majority hot electrons. This implies that most of the hot electrons 
crossing the DW will be majority hot electrons. However, once the 
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Figure 2 | Time resolved Kerr signal measured for different laser fluences 
on (a) the CoPt mulitlayers and (b) on the CoPd multilayers for various 
laser fluences (black 4.9 mj/cm^, red 4.2 mj/cm^, green 3.5 ml/cm^, blue 
2.8 mj/cm^, cyan 2.1 mj/cm^, magenta 1.4 mj/cm^, yellow 0.7 mj/cm^, 
dark yeUow 0.35 mj/cm^). The colored arrows schematically represent the 
increase in laser fluence with blue and red corresponding to low and high 
fluence, respectively, (c) Evolution of the demagnetization time as a 
function of the quenching of the Kerr signal for the CoPd (black squares) 
and CoPt (red circles) multilayers. 



DW is passed, majority hot electrons suddenly become minority hot 
electrons in the neighboring domain. Since they have much smaller 
lifetime and velocity, they will induce a minority spin accumulation 
in an area close to the DW, leading to a local reduction of the 
magnetization. 
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Figure 3 | (a) Time dependance of the Kerr signal measured on the CoPt and CoPd multilayers measured for different applied fields for a laser fluence of 
4.9 mj/cm^. (b) Dependance of the demagnetization time as a function of the initial magnetic state represented by the ratio of the Kerr signal 0^ obtained 
for an applied field HqH over the Kerr signal measured at saturation for the CoPt (red fiUed circles) and CoPd (open diamond) multilayers for a laser 
fluence of 4.9 mj/cm^. The red (black) dashed line corresponds to the estimated demagnetization time (see data analysis in the Methods section). 



In their paper, Pfau and collaborators evaluated the spin accu- 
mulation at the domain wall boundary using the superdiffisuve 
model'^'^'^. It was shown to result in a broadening of the magnetiza- 
tion profile around the domain wall center with a spatial extent of 
20 nm resulting in the modification of the DW scattering form factor 
and thus leading to a shift of the maximum scattering intensity. 
Therefore, the demagnetization signal and a possible alteration of 
the demagnetization time due to the hot electron spin transport 
between neighboring domains should only come from 8 = 20 nm 
around the DW. 

The dynamics due to hot electron spin transport should alter the 
demagnetization curve shape and result in an increase of the demag- 
netization at the peak position. This extra demagnetization has been 
shown to vanish within the electron-phonon relaxation time (about 
1 ps)"". However, no modification of the dynamics can be observed 
in the present study (Fig. 3a). Taking into account the noise level in 
our measurements, we can estimate the maximum proportion of the 
signal due to superdiffusive transport. The noise level can be esti- 
mated to be 3.5% and 5% for the CoPt and CoPd multilayers 



measurements, respectively. Using equation 4 (the supporting ana- 
lysis to equation 4 can be found in the "Methods" section), the 
demagnetization signal due to superdiffusive spin currents are smal- 
ler than 33% and 25% of the total demagnetization signal within the 
20 nm around the DW for CoPd and CoPt multilayers, respectively. 

In the same way, we can estimate the contribution of hot electron 
spin transfer on the total demagnetization time. With the domain 
sizes in our CoPd and CoPt multilayers being 132 and 145 nm 
respectively, and assuming that hot electron spin transport result 
in a constant demagnetization time of z^c = 100 fs in both systems^", 
the reduction of the effective demagnetization time between the 
saturated and the multidomain states would be 19 fs and 18 fs for 
CoPt and CoPd respectively, as represented by the dashed lines in 
Figure 3b (see data analysis in the Methods section for more details). 
In both cases, the expected decrease of the demagnetization time due 
to hot electrons transport is much larger than the error bars. We can 
thus conclude that the contribution of hot electrons spin transport on 
the demagnetization is negligible in our systems. As a result, the time 
resolved signal we measure by time resolved magneto-optical Kerr 
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effect is mainly reflecting demagnetization due to spin flip pro- 
cesses. Similar arguments can be used to analyze magnetic scattering 
experiments if only the amplitude of the signal is taken into 
account'". Even though, in the study of Vodungbo et al. the samples 
show a magnetic structure with domain sizes of 80 nm, using sim- 
ilar arguments as above, an increase of the demagnetization time 
with the laser fluence is still expected, contrary to what has been 
observed. Furthermore, our results clearly contrast with recent 
theoretical works supporting spin transport as the only relevant 
process in ultrafast demagnetization" '". 

Our conclusions also contrast with experiments realized on mag- 
netic spin-valve like multilayers for which a clear evidence of an 
ultrafast spin transport resulting in a modification of the magnetiza- 
tion dynamics was observed"' "'^'''^'. However, in that case, the dif- 
ferent magnetic layers are usually only few nm thick comparable with 
the inelastic mean free path of hot electrons'^. The spin accumulation 
due to hot electrons spin transport takes place over the full thickness 
of the layers inducing a demagnetization signal of the same order as 
the one generated by spin-flip scattering processes. 

Finally, we would like to stress that even though many theoretical 
and experimental works have focused on hot electron spin lifetime in 
magnetic materials''"", the continuous variation of the local spin 
orientation in the D W induces a rotation of the exchange field which 
might significantly influence the spin lifetime. For instance, it was 
recently shown that the spin diffusion of conduction electrons by a 
magnetic DW strongly depends on the DW structure, and more 
precisely on the magnetization gradient present over distances com- 
parable to characteristic lengths such as the spin diffusion length". 
However, while they could be of prime importance in determining 
the spin accumulation length around the magnetic DW, the effects of 
a magnetization gradient on the diffusion of hot electron spins 
remain an unexplored area. 

To conclude, we studied the demagnetization dynamics in CoPt 
and CoPd multilayers as a function of their magnetic configurations. 
When the samples are magnetically saturated, their demagnetization 
dynamics is similar to what is usually observed in 3d metals, namely, 
an increase of the demagnetization time is witnessed when the loss of 
magnetization increases. By measuring the demagnetization dyna- 
mics as a function of the domain structure, the extracted constant 
demagnetization time measured for different magnetic configura- 
tions reflects the absence of a modification of the demagnetization 
dynamics by hot electron spin transport in our samples. The analysis 
of the demagnetization curves shows that only a small part of the 
signal can be attributed to inter-domains hot electron spin transfer. 
Altogether, we showed that the role of inter-domains hot electrons 
spin transport on the demagnetization time is negligible. Reducing 
the domain size down to a length comparable to the spin accumula- 
tion length could potentially lead to a modification of the demagnet- 
ization dynamics due to inter-domain hot electron spin transport. 

Methods 

Sample preparation. The samples used in this study consist of two different 
multilayers SiN 30/ Al 10/[Co 0.4/Pd 0.6] x2o/Pd 2 and SiOx/Ta 5/Pt 5/[Co 0.5/Pt 
1-0] xis/Pt 5 (thicknesses in nm) grown at room temperature by dc magnetron 
sputtering. 

Experimental methods. The static polar Kerr loops presented in Fig. 1 were 
measured at a wavelength of 543 nm. The pump-probe measurements were carried 
out using a Femtolaser XL 300 pulsed laser with a wavelength of 800 nm and a pulse 
duration of about 50 ± 5 fs at a repetition rate of 5 MHz. The probe beam is 
frequency doubled (400 nm) in order to minimize optical artifacts such as dichroic 
bleaching'' \ The pump is much stronger than the probe with an intensity ratio of at 
least 20 for the lowest pump fluence. The pump laser is focused down to a spot size of 
about —30 fim while the probe beam is focused down to about —20 |im. Both beams 
are impinging on the sample almost at normal incidence with an angle of about —5 
degrees between them. In this configuration, the measured Kerr rotation is 
proportional to the out-of-plane component of the magnetization. The pump fluence 
is in the range of 0.5 to 4.9 mj/cm*^. 



Data analysis. The spin dynamics submitted to a laser pulse excitation can be 
described by a phenomenological thermodynamic model, i.e. the so-called three- 
temperatures modeV'^^. This model supposes that the system is split in three baths: the 
electron, lattice and spin baths. The laser photons are directly absorbed by the 
electronic system instantaneously creating hot electrons. Once thermalization is 
achieved by electron-electron interactions, each bath can be described by its own 
thermodynamic temperature. Then, the relaxation takes place by energy transfer 
between the different baths to balance their temperatures. The set of differential 
equations describing the heat flow between the three sub-systems are solved 
neglecting the spin specific heat, assuming that the measurements are done in the low 
laser fluence limit and that the electron temperature rises instantaneously upon laser 
excitation: 



= [( 



(A2T£-AiiM)e 



T£(Ai-A2)e"^ 



(1) 



)0(f)]®r(() 



where Y{t) is the Gaussian laser pulse, &{t) is the step function and b{t) is the Dirac 

delta function. The constant Ai represents the value of — after equilibrium 

between electrons, spins and lattice is restored. The constant A2 is proportional to the 
initial electron temperature rise, t^^ and T£ are the demagnetization time and the 
relaxation time, respectively. 

Assuming that the spin flip processes and hot electron spin transport contributions 
to the demagnetization time are independent we can calculate the effective demag- 
netization time as: 



11 1 



(2) 



where T^^and isc are the demagnetization times due to spin-flip processes and hot 
electron spin transport, respectively. In such case and for t^"^' — 200 fs and 
^coPd = fs as measured in the saturated state (Fig. 3b) and x.^ = 100 fs", the total 
demagnetization time Tm is given by: 



{d — d)*1s} + &* Xeff 

' d 



(3) 



with d being the domain size and 6 — 20 nm as the spin accumulation spatial extent 
around the domain wall. We can therefore estimate the demagnetization time for 



CoPt, 



aSl fs, and CoPd, rV; 



- 1 62 fs, multilayers in the multidomain state. 



Thiscorresponds to a reduction of 19 and 18 fsfor CoPt and CoPd compared with the 
demagnetization time measured in the saturated state. 

Similarly, the total demagnetization due to both spin flip processes and spin 
dependent transport can be expressed as: 



{d-3) * 3Msf-\-3 * (3M,f + 3M,c) 



(4) 



where tW^^and 3Msc are the demagnetization amplitudes due to spin flip processes 
and hot electrons spin depedent transport, respectively. 
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